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A  STUDY  OF  THE  MOTION  AND  STABILITY 


OF  TORPEDOES  IN  3  DEGREES  OF  FREEDOM 


1 .  INTRODUCTION 


In  this  paper,  the  dynamical  elements  which  are  necessary  for  the 
analytical  development  of  motion  and  stability  of  a  torpedo,  or  slender  body 
are  presented. 


Prom  classical  mechanics  it  is  well  known  that  the  translational 
motion  of  a  body  is  described  by  the  equation  of  linear  momentum  while  the 
rotational  motion  is  given  by  the  equation  of  angular  momentum.  In  general 
the  translational  and  rotational  motion  of  a  torpedo,  or  slender  body  are 
coupled  and  need  to  be  considered  simultaneously.  The  general  treatment  of 
this  problem  is  complicated  since  the  motion  is  described  by  six  degrees  of 
freedom.  It  is  possible  however  to  obtain  criteria  for  stability  for 
example,  in  the  horizontal  plane  by  considering  just  three  degrees  of  freedom, 
that  is,  by  considering  the  motion  of  the  centre  of  gravity,  parallel  and 
orthogonal  to  the  longitudinal  axis  of  the  torpedo  and  rotational  motion  about 
a  horizontal  transverse  axis  through  the  centre  of  gravity.  A  similar 
procedure  can  also  be  used  to  derive  the  equations  of  motion  in  the  vertical 
plane.  These  equations  are  all  highly  non-linear.  However,  on  linearizing 
and  decoupling  the  motion,  a  set  of  ordinary  differential  equations  for  the 
pitch  and  yaw  angles  are  found.  The  Laplace  transform  of  these  differential 
equations  yield  criteria  for  stability  in  the  vertical  and  horizontal  planes. 


In  the  case  of  a  torpedo  moving  through  fluid  at  rest  at  infinity, 
the  linearized  equations  can  also  yield  details  of  the  trajectory.  Analysis 
of  the  torpedo  trajectory  will  be  discussed  elsewhere  [1], 


2.  EQUATIONS  OF  MOTION  IN  THE  VERTICAL  PLANE 


The  equations  of  motion  discussed  in  this  section  describe  the 
dynamical  behaviour  of  a  torpedo  in  the  vertical  plane,  moving  with 


rectilinear  motion  through  a  stationary  ideal  fluid  of  infinite  extent.  The 
question  of  dynamic  stability  and  the  derivation  of  the  equations  of  motion 
will  be  addressed  under  the  following  assumptions: 


(a)  the  thrust,  velocity  and  weight  of  the  torpedo  are  constant, 
with  the  thrust  tangential  to  the  trajectory, 

(b)  the  shape  of  the  torpedo  is  taken  to  be  a  prolate  spheroid  of 
finesse  ratio  l/d  (£  >  d),  where  l  is  the  length  of  the 
torpedo  and  d  its  diameter;  and 

(c)  cavitation  and  free  surface  effects  are  assumed  to  be 
negligible. 


Hie  reason  for  assuming  the  torpedo  is  a  prolate  spheroid  is  to 
enable  the  mass  accession  (or  added  mass)  due  to  the  motion  of  the  torpedo 
through  the  fluid  to  be  incorporated  into  the  dynamical  equations  of  motion. 
Moreover  when  a  body  accelerates  in  a  fluid,  then  not  only  does  the  torpedo 
accelerate  but  the  fluid  around  the  torpedo  generally  also  moves.  This  means 
that  the  torpedo  has  to  do  work  to  increase  the  kinetic  energy  of  both  the 
torpedo  and  the  fluid. 

Hie  nomenclature  to  be  used  throughout  this  paper  in  writing 
equations  of  motion  is  that  the  translational  velocities  will  be  denoted  by 
components  (u,v,w)  and  the  rotational  velocities  by  (f,&,$).  (’)  denotes 

differentiation  with  respect  to  time. 


A  torpedo  in  rectilinear  flight  is  subject  to  a  number  of 
hydrodynamic  and  hydrostatic  forces  as  a  consequence  of  the  water  pressure 
acting  externally  on  the  torpedo.  A  force-moment  diagram  is  shown  in  Figure 
2.1,  for  a  torpedo  on  a  steady  course  at  set  depth  in  the  vertical  plane.  In 
this  diagram  T  is  the  thrust  generated  by  the  propulsion  unit  in  newtons  (N) 
and  is  directed  along  the  longitudinal  axis  with  no  resultant  torque  or 
misalignment.  L  and  D  are  the  resultant  lift  and  drag  forces  in  newtons 
respectively,  acting  parallel  and  orthogonal  to  the  velocity  V  (m/s)  in  the 
(x,z)  plane.  B  is  the  buoyancy  (N)  and  represents  the  nett  force  of  the 
hydrostatic  pressures.  W  is  the  weight  (N),  and  MQ  the  measured  moment 
(N.m).  N&  represents  the  damping  moment  which  is  proportional  to  the  angular 
velocity  and  F&  is  the  transverse  force  associated  with  the  damping  moment  in 
the  vertical  plane,  a  is  the  orientation  angle  in  pitch,  3  is  the  attack 
angle  in  pitch  and  3  is  the  elevator  angle  in  the  vertical  plane. 


In  order  to  find  the  equations  of  motion  for  a  torpedo  moving 
through  an  incompressible  fluid,  we  firstly  derive  the  kinetic  energy  of  the 
fluid.  When  the  velocity  _u  of  the  fluid  is  irrotational,  we  may  introduce  a 
velocity  potential  <p  satisfying: 

u  *  , 


2 


0 

FIGURE  2-1  Hydrodynamic  and  hydrostatic  forces  acting  or  torpedo 
(i)  x,yfz  Body  Axes 

Hi)  x'.v'.z'  Absolute  Reference  Frame 


and  because,  for  incompressible  flow,  (7.u  «  0)  it  follows  that 


2  * 

V  ♦ 
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* 

For  an  irrotational  flow  characterized  by  the  potential  $  ,  the  kinetic  energy 
Ep  in  the  fluid  domain  A  is  given  by  [2]: 


i 


-f/\  *"  <» 

A 


1,2,3) 


where 


denotes  differentiation  of  $  with  respect  to  u 


Equivalently: 


ep  -  -  if  ♦*V‘ds- 

s 


(2.1 ) 


where  jn  is  a  unit  normal  to  the  surface  S  of  a  slender  body  and  p  the  density 
of  the  fluid.  Let  the  translational  velocity  components  of  the  torpedo  or 
slender  body  be  given  by: 


(u1 , u2,u^)  •  (u,v,w) , 


and  the  rotational  motion  by: 


(u4,u5,u6)  -  (-f,A,$). 

For  each  component  of  the  motion  there  exists  a  corresponding  velocity 
potential  for  the  fluid.  That  potential,  per  unit  uv,  will  denoted  by  <t>v. 
Hence 


♦  -  uv<frv;  (v  -  1 , — ,6)  ,  (2.2) 


i . e. , 


♦,u  +  $2v  +  <|>3w  +  +  ij>5&  +  <t>6?). 


(2.3) 


Substitution  of  (2.2)  into  (2.1)  yields: 


F  ‘  pv  11  V 


where  the  cartesian  tensor  A  is  defined  by: 


V  *  -p 


>  /  $  3  <|>  dS , 
J  n  v 


and  satisfies  the  symmetry  condition: 


A[uv]  *  0  * 


We  call  the  elements  of  A^y  the  components  of  the  added-*mass  tensor  [2]. 

When  there  is  only  motion  in  the  vertical  plane  (i.e.,  the  (x,z) 
plane)  the  velocity  potential  is  spanned  by  just  three  degrees  of  freedom  and 
may  be  written  in  the  absence  of  torpedo  roll  (non-rolling  coordinate  system) 


+  <fr3w  +  <fr5&  . 


(2.5) 


It  follows  from  (2.4)  that  the  quadratic  form  of  the  kinetic  energy  of  the 
fluid  in  this  coordinate  system  is: 


EF  ”1/2{A11u2  +  A3;jw2  +  2A35w&  +  A5542}, 


(2.6) 


with  all  other  components  of  A  zero  when  $  is  given  by  (2.5).  This  results 
from  the  symmetry  of  the  torpedo  ie.,  replacing  u  by  -u  should  leave  energy 
invariant;  to  A13  >  A^5  «  0. 

The  kinetic  energy  of  the  torpedo  can  similarly  be  expressed  as  a 
quadratic  form,  i.e.. 


E  »  V?  B  U  u  ;  B.  ,  *  0, 
T  *  yv  y  v  (yvj 


where 
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In  (2.7),  Ijj  are  the  components  of  the  moment  of  inertia  tensor,  xQ,  yQ  and 
zQ  the  coordinates  of  the  centre  of  mass  and  H  the  mass  of  the  torpedo.  For 
motion  in  the  vertical  plane  with  the  centre  of  gravity  of  the  torpedo  located 
directly  below  the  longitudinal  axis  (a  distance  zQ)  the  kinetic  energy  of  the 
torpedo  is  given  by: 


0  0  o 

Et  »y2{Mu  +  Mw  +  I224  +  2MZoU&^'  (2.8) 

and  hence  the  total  energy  by: 

E  M)u2  +  (A,,+  M)w2  +  (I  _  +  A__)&2  +  2A,_wd  +  2Mz  u4}  .  (2.9) 

11  33  2253  35  O 


If  and  _H  are  the  linear  momentum  and  total  angular  momentum  respectively 
then: 


G  ■  ej3  E  -  {(M  +  A. . )u  +  Mz  &,0,(M  +  A,,)w  +  A,_&}  ;  j-1,2,3,  (2.10) 

—  — -3  Uj  n  o  33  35 

and 


H 


E 


{0, (I  +  A  )&  +  Mz  u  +  A  _w,0}  ;  k=4,5,6  . 

ii  55  O  35 


(2.11) 


Furthermore  if  _F*  and  _L*  are  the  resultant  external  force  and  momentum 
respectively  of  both  the  torpedo  and  the  fluid  and  if  an  inertial  reference 
frame  is  fixed  to  the  torpedo  moving  with  velocity^  -  (u,o,w)  and  angular 
velocity  u>  *  (o,A,o)  then: 


F*  *  G  +  u)XG, 


(2.1 2a) 


and 


L*  «  H  +  0>XH  +  uXG  ,  (2.12b) 

where  X  is  the  vector  product.  More  explicitly: 

F*-{(M  +  A,,)u  +  MZ  O  +  (M  +  A-i-Owd  +  A__  d2,0, (M  +  A,,)w  +  A.,  a  -(M+A,,)ud  -  Mz  d2} 
—  II  O  35  33  35  1  1  O 

(2.13) 

and 

_L*  »  {0,(I22  +  h55)a  +  Mzqu  +  A jgW  -  uw(AJ3  -  A^)  -  A^ua  +  MzQwd,  0> . 

(2.  r.) 

Resolution  of  the  forces  in  Figure  2.1  acting  parallel  and 
orthogonal  to  the  direction  of  motion  yields: 

F*  -  {T  -  DcosB  +  (B-W)sina-Lain8,0,-Dsinf5-(B-«)cosa-LcosB  +  Fd} .  (2.15) 

Similarly,  the  resolution  of  torques  acting  about  the  centre  of  gravity 
yields: 

L*  -  {  0,M  -  Ad  +  Bx  cosa  -  Wz  sinarO}.  (2.16) 

—  o  o  o 

Equating  (2.13)  and  (2.15),  (2.14)  and  (2.16)  the  equations 
governing  the  notion  of  the  torpedo  in  the  (x, z)  plane  can  be  expressed  as: 

T  -  DcosB  +  (B-W)sina  -  LsinfJ  -  (M  +  A..)i  +  M*  a  +  (M  +  A,,)wd,  (2.17a) 

11  o  J 

2 

-DsinfJ  -  (B-W)cosa  -  LcosB  +  Fd  =*  (M  +  A,,)w  +  A  ,_a  -  (M  +  A,,)ud  -  Mz  d  , 

33  35  11  O 

(2.17b) 

M  -Ad  +  Bx  cosa  -  Wz  sina  *  ( I__+  A__)a  +  Mz  u  +  A,_w  -  uw(A. ,  -  A,,) 

O  O  O  22  55  0  35  11  33 

-A,_ud  +  Mz  wd,  (2.17c) 

35  O 

where  (u,o,w)  ■  (VcosfJ ,0,VsinB ) ,  and  V  is  defined  in  Figure  2.1. 
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LINEARIZED  EQUATIONS  OF  MOTION  IN  THE  VERTICAL  PLANE 


In  general,  analytic  solutions  to  (2.17)  are  not  possible,  however 
upon  linearizing  these  equations,  solutions  can  be  found.  To  do  this  we 
assume  that  the  hydrodynamic  forces  and  moments  are  approximately  linear  in 
attack  angle  and  turning  rate. 


The  forces  and  moments  appearing  in  (2.17)  are  both  functions 
of  3  and  6  which  can  be  taken  to  be  independent  and  small  (less  than  10°). 
Expanding  L  and  MQ  in  a  Maclaurin  series  and  neglecting  all  but  first  order 
terms  yields: 


S(3,5)  -  s(0,0)  +  s66  +  Sp3, 


where  S.  for  example  denotes  partial  differentiation  of  S  with  respect  to  6. 
If  the  torpedo  is  highly  symmetrical  we  can  assume  without  loss  of  generality 
that  S(0,0)  =»  0.  Furthermore  if  a  is  also  small  then  we  have: 


cosa  =  cosB  -  1 , 
sina=a,  sinS^B  . 


Incorporation  of  the  above  approximations  into  the  equations  of  motion  yield 
the  following  set  of  linear  equations 

(T  -  D)  +  (B  -  W)a  =  Mz  a  ,  (3.1a) 

o 

-(D  +  Lg)B  -  A^a  -  (B  -  W)  -  Lg6  +  F&  »  -mtv«S  +  ,  (3.1b) 

M.6  +  Mi  +  M  &  +  Bx  -  Wz  a  -  A,_  V§  *  Ea  ,  (3.1c) 

0  12  o  o  35 

where 


m 

T 

M  +  A^, 

(3.2a) 

mT  : * 

M  +  A,,, 

(3.2b) 

L 

33 

V" 

V  v2(»i r  a33>- 

(3.2c) 

M„  :  = 

-(A-VA,,) , 

(3. 2d) 

In  terms  of  the  coefficients  of  inertia  for  a  prolate  spheroid  [3,4] 
(3.2  a,b,c,e)  can  be  expressed  as: 


»T  “  M  +  k1Mf, 


“L  "  M  +  W 


mt  -  Me  -  v  Vki  -  V' 


£  -  i22  +  k’Mf' 


where  is  the  displaced  fluid  mass  and  k1 ,  k2  and  k'  the  coefficients  of 
inertia: 


k  -  o  (2-a  )  , 

i  o  o 


k  -  6(2-6)  , 

2oo 

k*  -  (a2-b2)  2(6  -a  )  { (a2-*2)  [2(a2-b2)  -  (a2-*2)  (6  -<x  )] 
o  o  o  o 

where  a  is  the  semi-major  axis  and  b  the  semi -minor  axis  of  the  prolate 
spheroid .  Furthermore : 


1  — e  /,  (1  +  e)  „  \ 

~r  \lmTT-r-^r  *  27- 

© 

1-e2  /  e  i,  ,_<1  +  e)\ 

T5"  (777 ' /2 


where  e  is  the  eccentricity. 

If  p  is  the  density  of  the  fluid,  and  A,  l  and  V  the  cross-sectional 
area,  the  length  and  velocity  of  the  torpedo  respectively,  then  (3.1b)  and 
(3.1c)  can  be  expressed  in  terms  of  the  hydrodynamic  coefficients  which  are 
defined  below  on  dividing  (3.1b)  by: 


A1 :-  V2  PAV  , 


(3.3a) 


and  (3.1c)  by 


together  with  a  change  in  the  independent  time  variable  t  ie. ,  the  time 
variable  t  is  changed  to  the  dimensionless  arc  length  parameter  C , proportional 
to  t  and  defined  by: 


where  tc  is  called  the  characteristic  time.  On  substitution  of  (3.3a)  and 
(3.3b)  into  (^.1b)  and  (3.1c)  respectively  it  becomes  apparent  that  tc  is  best 
chosen  as  IV~  ,  where  l  is  the  length  of  the  torpedo.  This  yields: 


I 

and 


-  A*  a* '  *  (C  +  m*)a' 
35  L 


mB*  +  Ct  )3  +  C  6 

T  D  Le  Ls 


B  >  0,  (3.4a) 


£*«••  -C^a'  .y'  -  V  -CMS 


**  ** 

B  +  W  az  -  0, 
o 


(3.4b) 


where  ()'  denotes  differention  with  respect  to  5.  In  the  vertical  plane  the 
hydrodynamic  coefficients,  corrected  for  the  effects  of  entrained  water  are: 


-VT'  • 

-  fO^pav*)-1 

-  MiA"1  , 

:-«2(1/2PAVA2) 

-  L(JA1  , 

■  vl1  - 


-1 


is  the  drag  coefficient, 
the  lift  coefficient, 
the  transverse  force  coefficient, 
the  moment  coefficient 
the  damping  moment  coefficient 
the  elevator  lift  coefficient, 
the  elevator  moment  coefficient. 


Other  terms  appearing  in  (3.4)  are  defined  as  follows: 


A  :-  A  (VfcpAf2)”1;  m  =  M  0/2  pfJL )  ~ ' ;  m"=  m  0/2  P**>  _1  ; 


-1 


-1 


35  35 


T  T 


*  — 1  *  ,  3  -1  **  -1 

B  :»  (B-W)A1  ;  I  :=  E C/2  PAi  )  }  B  :*  X  BA„ 


O  2 


**  -1 
W  :*  WA2  . 


4.  EQUATIONS  OF  MOTION  IN  THE  HORIZONTAL  PLANE 


The  equations  describing  the  motion  in  the  horizontal  plane  can  be 
derived  following  a  similar  procedure  to  the  one  given  for  deriving  the 
equations  in  the  vertical  plane.  A  number  of  minor  changes  are  all  that  are 
necessary  so  not  all  the  details  are  given  here.  For  example,  in  the 
horizontal  plane  the  gravity  terms,  weight  and  buoyacy  can  be  equated  to  zero 
in  both  the  force  and  moment  equation. 


Take  S  to  be  the  rudder  angle,  $  the  orientation  angle  in  the 
horizontal  plane  (i.e.  yaw  angle)  and  Mf  the  angle  of  attack  in  yaw. 
Furthermore  u.  ■>  (0,0,$)  and  ji  =  (u,v,o).  With  no  displacement  of  the  centre 
of  gravity  in  the  (x,y)  plane,  the  total  kinetic  energy  of  both  the  fluid  and 
torpedo  can  be  expressed  as: 


E  -V2{(M  +  An)u2  +  (M  +  A22)v2  +  (A66  +  133)$2  +  2A26v$}.  (4.1) 


From  (4.1)  the  linear  momentum  and  total  angular  momentum  have  the  form: 


£»  {(M  +  An)u,  (M  +  A22)v  +  A26$,  0}, 

and 

JL  *  {o,0,(Agg+  *33^  +  A26vf ' 


respectively.  From  (2.12)  it  further  follows  that  the  external  force  and 
momentum  L*  are  given  by: 


F.*  -  { (M  +  An)u  -(M  +  A22)v$  -  A26$2,  (M  +  A22)v  +  A2g$  +  (M  +  An)u$,  0}  ,  (4.2a 


and 
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_L*  -  {0,0,  (Agg  +  I33)  4>  +  A26(v  +  i >u)  +  uv(A22  -  An)}, 


respectively.  Expanding  the  left  side  of  (4.2)  in  a  Maclaurin  series, 
neglecting  all  but  first  order  terms  gives: 


F*  -  (F1 (u, $  ,6 ) ,  F2 (v,| ,5 ) ,0)  = 


(F^u  +  fU  +  F§5,  F2v  +  F?<fr  +  F26,0),  (4.3a 


and 


L*  -  (0,0,LJ(v,$,6) )  *  (0,0, v  +  L?i  +  L|5)  . 


Equating  (4.2)  and  (4.3)  yields: 


F^u  +  fU  +  F$S  “  mjU  -  ny4  -  A26<fr2, 

♦ 


F^v  +  F?$  +  Fj6  »  nVpV  +  A26^  +  mLu$ 


L^v  +  lH  +  L$$  *  0$  +  A2g{'fr  +  +  uv(A22  -  An>  f 


where 


M  +  A. 


“l  "  "  T  "11  ' 

■  M  +  A. 


22  ' 


0 


I33  +  A66  ' 


and 


(u,v,o)  «  (Vcos'f,-  Vsin^^),  V  *  constant. 


(4.3b. 


(4.4a) 

(4.4b) 

(4.4c' 


5.  LINEARIZED  EQUATIONS  OF  MOTION  -  HORIZONTAL  PLANE 

Following  the  method  of  Section  3,  equations  (4.4b)  and  (4.4c)  can 
be  expressed  as: 


1 


+  (<V  'V*'  +  V 


CM  ?  +  CM. 


+  A26HM  +  \* 


A26^’  +  V' 


-  0  $"  »  0. 


(5.1a 


(5.1b) 


In  the  (x,y)  plane  the  hydrodynamic  coefficients  corrected  for  the  effects  of 
entrained  water  are  defined  by: 


Fv  V 

%  ' 

Cp  :«  F?  ^pAVt)-1  , 
t  ♦ 

„  _2  -1 

CLfi  S“  F6  A1 

„  T  3.-1 

CM6  !“  V2  ' 

Cw  5*  L.^PAVt2)-1  , 
Md  2 

CM  I'/ZiPAV2*)'1  , 


the  lift  coefficient. 


transverse  force  coefficient. 


rudder  lift  coefficient, 


rudder  moment  coefficient. 


deunping  moment  coefficient. 


moment  coefficient. 


where  L.,  :-  VL2  -  V2(A22-A11),  L2:«  L?  -  A2gV,  A2g:«  A^t1^  PAJt2)-1 

$ 

*  ji  3  — 1  *  * 

and  0  0(72  PAf.  )  .  Hie  coefficients  mT  and  m^  have  the  same  algebraic 

form  as  given  in  Section  3. 


6.  STABILITY  CRITERIA  FOR  MOTION  IN  THE  HORIZONTAL  AND 


VERTICAL  PLANES 


Hie  problem  of  determining  stability  of  a  linear  system  is  one  of 
finding  roots  of  the  characteristic  equation,  found  on  taking  the  Laplace 
transform  (L-transform)  of  the  equation  governing  the  motion.  In  general  if 
the  characteristic  equation  of  the  transformed  equation  is  given  by  [51: 


o  sn  +  a. s11”1  +  ...  +  a  , s  +  a  , 
o  1  n-1  n 


(6.1 ) 
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then  the  criteria  for  stability  are  that  the  coefficients  of  the 
characteristic  function  (6.1)  have  roots  that  lie  in  the  left  half  of  the  s- 
plane.  This  is  both  a  necessary  and  sufficient  condition  and  is  equivalent 
to  saying  all  the  determinants: 


D1-V  D2 


must  be  positive,  i.e.. 


>  0,  >  0,  Dj  >  0,  ...  • 


If  any  coefficient  of  any  power  of  s  is  either  negative  or  zero,  then  there  is 
at  least  one  root  in  the  right  half  of  the  s-plane  or  on  the  imaginary  axis  of 
the  s-plane  in  which  case  the  system  is  unstable.  A  more  complete 
description  of  the  motion  is  unnecessary  since  all  that  is  required  is  a 
knowledge  of  whether  or  not  the  system  is  stable,  i.e.,  whether  a  response  to 
a  perturbation  remains  bounded  or  grows  as  t  (or  ()  +  •. 


(a)  Conditions  for  stable  motion  in  the  horizontal  plane 


On  differentiation  of  (5.1)  it  is  not  difficult  to  show  that  the 
decoupled  equation  for  yaw  angle  $  in  the  €  domain  can  be  expressed  as: 


+  a,$"  +  a_$'  +  a, 6  +  a  5  -  0, 

4  3  2  10 


(6.2) 


where 


*  *  *2 
°4  *"  "T  °  '  A26  ' 


a3  !”  ~®T  CMd  +  A26(CFt  “  "l5  “  CLI?0  “  °MA26  ' 


a.  -(m_  -  C_  )C..  +  C.  C..  , 


(6.3) 


i 


'OTJ','7  W7^v.V  V-'  ;.'A  ..■%  ■>  .->-f.%  v; 


S 


Without  loss  of  generality  we  assume  that  the  desired  heading  is  $  *  0,  then 
(6.2)  represents  the  error  in  heading.  For  stability  it  is  necessary 
that  ♦♦0. 

The  L-transform  of  (6.2)  yields: 

3  2  2 

♦  («V  +  a,s  +  a.s)  +  S(a,s  +  a  )  «  (s  ♦  a(  +  s(a.<J>'  +  a,*  )  +  a.**' 

4  3  2  1  O  04  4030  40 


♦v; 


+  a  6  +  a  6  ) 

2  o  1  o 


(6.4) 


where  s  is  the  transform  variable  and  4  ,  6',  4 ' '  6  ,  initial  conditions. 

o  o  o  o 

The  stability  of  the  torpedo  in  the  horizontal  plane,  can  be 
inferred  provided  the  rudder  response  6  is  known.  For  example  if  c.|  is  a 
real  constant  ie.,  6  »  c^6,  then  (6.4)  becomes 


♦  (s) 


s  ♦  a.  +  s(a  ♦'  +  o,6  )  +  a. 6''  +  a  6*  +  a  6  +  a  6 

04  40  30  40  30  20  IQ 

3  2 

as  -fas  4*  (a  +  ca)s+ac 
4  3  2  1  1 7  o  1 


P(s) 
Q(s)  * 


For  stability  the  degree  of  the  polynomial  Q(s)  must  be  greater  than  P(s)  and 
the  roots  of  Q(s)  in  the  complex  a -plane  should  lie  in  the  left  half. 
Furthermore  the  determinants: 


3  ; 


should  all  be  positive. 


°2  - 


a  c. 
o  1 


°2  +  Clai 


°3  ‘  D1°2  ' 


If  the  torpedo  has  neutral  rudder  response,  that  is  S»0  and  further 
if  the  ofigin  of  the  coordinates  is  at  the  centre  of  the  prolate  spheroid, 
(i.e.,  *26”  0  ••  *•  C1  *  then  the  above  determinants  reduce  to: 


o3>0  and  a2>0  . 


(6.5) 


For  small  angles  of  attack  the  stability  criteria  for  a  torpedo  without 
controls  can  be  expressed  in  terms  of  the  hydrodynamic  coefficients  as: 


and 


V 


*  < 


(6.6a) 


at 
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CC  -  (m_  -  C  )  C  >  0  , 

m,l,b  vl  f.  m  ' 


The  general  solution  o£  (6.2)  for  neutral  rudder  response  and  «  0  can  be 
expressed  as: 

a  C  a  C 

<p(i)  -  Ae  +  Be  “  , 


where 


a±  =  - 


2  Vo 

a  *  (a  -  4a  a  )'* 
3  *  {  3  W2  4J 


are  the  roots  of  the  auxiliary  equation  of  (6.2).  For  dynamic  stability 
<fr(C)+0,  which  means  the  roots  a  must  have  negative  real  parts,  as  well 
as  a^<4a^a  .  This  will  be  the  Case  only  if  (6.5)  is  satisfied,  since  is 
always  positive. 

A  second  case  of  interest  refers  to  the  motion  of  a  torpedo  with 
zero  attack  angle  in  yaw  (i|i*0)  and  no  controls.  In  such  cases  the  above 
stability  analysis  yields  CM  >  0  for  the  transverse  damping  moment 
coefficient.  ™d 


(b)  Condition  for  stable  motion  in  the  vertical  plane 


Following  a  method  similar  to  the  one  described  for  stability  in  the 
vertical  plane,  on  decoupling  (3.4)  the  equation  for  the  pitch  angle  a  is: 


B,a"*  +  B_a"  +  B  a1  +  3„a  +  B„6  -*-3,6+3  ■  0  ,  (6.7) 

6  5  4  3  2  1  0 


where 


*2  *  * 

06!“  A35  +  ®TE  ' 


A35<CF.+  ®L)  "  mTCM.  “  A35CM  +  1  (CD  +  CLft}  ' 
t  a  p 

*«-  "i  -  CM(CF  *  V  -  CM  'CD  *  CL,  >  ’ 

t  a  p 


3.:-  W  z  (C_  +  C.  ) 


•  •  •  ■ 


*  * 

8  •*  A _ C  —  m  C  # 

2  35  Lj  T 


ft  :*  -C  C  -C  (C  +  C  ), 
1  ML5  Lg ' 


V*  B*CM  -  B*#(CD  +  CLe}  • 


As  in  case  (a),  stability  criteria  can  be  found  if  the  response  6  is  known. 
For  example  if  6  «  C2a'  t*Je  L-transform  of  (6.7)  yields: 


s3o  0  +  s2(a  6  +  a'0  )  +  s(a"6  +  a  *  0  +  a  0  +  c.a  0_)-0 

O  O _ O  5 _ O  O _ O  6  0  3 _ o  4 _ 2  O  2  O 

s{06s3  +  05s2  +  s(04  +  c2^2>  +  ^3  +  ^  i  °2^ 


and  a  necessary  condition  for  stability  is  that  all  the  determinants 


°i  -  e5  , 


d2- 


8  +  c  6 

3  2P1 


°3  '  D1D2  ' 


S4*  °282 


be  positive. 


7.  DISCUSSION 


The  model  described  above  considers  some  basic  aspects  of  torpedo 
motion  and  stability  in  both  the  vertical  and  horizontal  planes.  Hie 
question  of  accessibility  of  the  stability  criteria  depends  on  a  knowledge  of 
the  hydrodynamic  coefficients  for  a  particular  torpedo  body.  In  the 
foregoing  discussion  no  explicit  studies  were  made  of  torpedo  stability  due  to 
lack  of  data  relating  to  hydrodynamic  coefficients  of  the  torpedo  to  which 
this  report  was  directed.  Therefore  this  study  should  be  considered  as  a 
first  step  only  in  developing  a  knowledge  of  motion  for  submerged  bodies. 
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